Preprint typeset in JHEP style - HYPER VERSION 



Imperial /TP/09/D J W/0 1 



Quantum kink and its excitations 



Arttu Rajantie, David J Weir 

Theoretical Physics, Blackett Laboratory, Imperial College London SW7 2AZ, UK 



E-mail: 



a.rajantieOimperial . ac .uk, lavid.weirOSOimperial . ac .uk 



Abstract: We show how detailed properties of a kink in quantum field theory can be 
extracted from field correlation functions. This makes it possible to study quantum kinks 
in a fully non-perturbative way using Monte Carlo simulations. We demonstrate this by 
calculating the kink mass as well as the spectrum and approximate wave functions of its 
excitations. This way of measuring the kink mass has clear advantages over the existing 
approaches based on creation and annihilation operators or the kink free energy. Our 
methods are straightforward to generalise to more realistic theories and other defect types. 
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1. Introduction 

Kinks, domain walls and other topological defects are of interest across physics. The 
existence of magnetic monopoles, for example, is a general prediction of grand unified 
theories . Cosmic strings Q may give an observable contribution to the cosmic microwave 
background radiation ^ or gravitational wave background [Q]. In string theory, D-branes 
are analogous to topological defects. 

There has been widespread interest in the excitation spectra of topological defects, 
although mostly as a calculational tool for finding perturbative corrections to the mass 
of the defect [§, §, IC, 11 1- However, excitations themselves are of physical in- 



terest 1 12, 13, 14 1 . Bound states of point-like defects correspond to new particle species, 
and in the case of extended defects they would be particles propagating on the brane. 
For instance, in braneworld models all the standard model particles can be thought of 
as states localised on a domain wall [^]. Finding stable particles localised on defects in 
fully non-perturbative approaches would therefore correspond to completely new particle 
species; indeed, for strongly-coupled theories there is no a priori reason to expect that the 
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perturbative spectrum is accurate. The bound states of cosmic strings have been previ- 



ously studied as a potential way of stabilising string loops |1£, 17, 18|. In some sense, 
reconstructing the excitation spectrum is the opposite approach to that taken classically 
in Vachaspati's paper on the reconstruction of a field theory [p^. 

Domain walls, vortices and strings are also of interest in condensed matter systems. 



Charge carriers lying on domain walls have been posited as a conduction mechanism |20|. 



Interfacial phenomena such as the wetting of a surface can be viewed as the motion of a 



domain wall in a curved spacetime [2^, 22 1. Kinks are also studied in integrable systems 
and in non-Hermitian field theories [p^ ]. 

In this paper we investigate field correlation functions in the presence of a kink, and in 
this sense our approach bears some similarity with Refs. [12, ^ 0]. We show that the cor- 



relators can be used to measure non-perturbatively the mass of the kink and its excitation 
spectrum. Furthermore, we obtain approximate wavefunctions for the excitations based on 
Monte Carlo data. Previously these have only been calculated at the linear level, although 
there has been some recent progress in accounting for interactions using renormalisation 
group techniques p^]. 

Previous work has measured the kink mass with 'Kadanoff's non-local operator' p6|| ; 
but more generally, the mass of topological defects has usually been calculated on the 



lattice by finding the free energy of the defect |27, 3C]. This is taken as the 

difference in free energy between the system with the defect present, and the same system 
but without the defect. The defect is often created using twisted boundary conditions, so 



one is effectively measuring the response to that twist [31|. Since the partition function 
cannot be determined by Monte Carlo simulations, one must measure derivatives of the 
mass and then perform finite differencing from a point where the mass is known. Typically, 
this is at a phase transition which is where the errors in Monte Carlo simulations are at 
their largest. 

The method of calculating the kink mass we propose requires study only of a topo- 
logically nontrivial sector, and only at one parameter choice. This gives a check on errors 
which may be present in the 'twist' method. It may also serve to reduce the amount of 
calculation required since it only requires Monte Carlo simulations to be carried out for 
the desired parameters. 

Our analytical discussion of correlation functions in the presence of topological defects 
takes place in Section |3[ Section 3.1 outlines methods which rely on the spectral expansion. 



We describe how to calculate the mass in Section 3.2, while in Section 3.3 we determine 
the particle spectrum of X(j)^ theory and the approximate wavefunctions in the presence 
of a kink. Monte Carlo simulations and other numerical work based on these ideas are 
discussed in Section ^. 

We believe that these techniques can be generalised to more sophisticated defects such 
as cosmic strings and monopoles; the extension to domain walls is trivial. Eventually, 
similar techniques may be used to study the defects present in supersymmetric theories 
such as the Wess-Zumino Model or SQCD and hence as a test of lattice supersymmetry p2|. 
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1.1 The X(j)^ model 

Kinks, or their higher dimensional counterparts domain walls, can occur where the vacuum 
manifold after spontaneous symmetry breaking is not simply connected, and the field can 
take on a different vacuum expectation value in different parts of the space. Domain walls 
could have formed in the early universe due to random fluctuations leading to different 
patches of space occupying different, disconnected vacua. 

In this paper, we focus on kinks in the 1+1-dimensional A(/>^ model, with the Lagrangian 



2! 



4! 



(1.1) 



with > in the classical broken phase. The vacuum is then 



±my^6/X. In (1 + 1) 



dimensions, the dimensionless parameter that appears in perturbative results will be Xjrr? . 

A kink configuration interpolates between the two vacua. We can create a topological 
kink by requiring that — > (/>o as x — > +oo and ^ —(j)Q as x ^ — oo, with one choice of 
sign being termed a 'kink' and the other an 'antikink'. The kink solution is obtained via 
the Bogomol'nyi equation for the kink H, 



The static classical kink solution (j)]^ is a minimum-energy configuration satisfying (|L 
that interpolates between the two vacua in the broken phase. In the frame of the kink. 



(1.2) 



4>]i{x, t) = m\j — tanh ( 



/ mx 

V7f 



The classical mass M^i is, from integrating the energy density. 



Mel = 4^/2^. 



(1.3) 



(1.4) 



1.2 Bound and scattering states 



We wish to consider the excitation spectrum in the presence of a kink. Expanding about 
the classical kink solution, (j) = + (j) and eliminating terms using the equation of motion 
for the classical kink yields [Ea] 



C 



3m 



— tanh — — 
2 \V2 



m 

T 



+ c 



int) 



leading to the form of the action, after integration by parts H 



(1.5) 



S 



cPx (j){x) 



-d^ + — 3m^ tanh^ 



mx 



(x) + 



int 



(1.6) 



where ^int are terms of order 0'^ and above. The Green's function G is given by 



m 



3m^ tanh^ 



mx 

7^ 



G{xi,ti;x2,t2) = d{xi - X2)d{ti - ^2)- 



(1.7) 
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We can construct this Green's function out of the eigenfunctions for the corresponding 
Schrodinger equation. Time-translation symmetry is not broken by the kink, so a suit- 
able complete set in the time direction is exp(i£'„i). On a timeslice, we must solve the 
Schrodinger equation with a Poschl- Teller potential: 



mx 



+ 3m tanh — — 
V2 



m 



with the eigenfunctions 



cosh 



mx 
mx 



. , / mx 
^ , smh — — 
V2J \V2 



/ iqmx\ 



3 tanh^ 



tpo = 

tpi = cosh~ 
ipg = exp 

and corresponding energy eigenvalues 

El 
E„ 



mx 

7f 



I — q — 3iq tanh 
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7f 
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(1.8) 

(1.9) 
(1.10) 
(1.11) 

(1.12) 
(1.13) 

(1.14) 



Eq corresponds to a Goldstone mode that shifts the kink in the x-direction; Ei corresponds 
to a massive localised 'bound state' and Eg are a continuum of scattering states labelled 
by — oo < g < oo. We impose antiperiodic boundary conditions on a box of length L with 
L — > oo. Then, the following expression for the allowed parameters q can be found, based 
on odd integer multiples of vr/L plus a phase shift due to the kink 



mL 



qn,a.nt-^ = (2n + 1)tT + S{qn,ant); 



n = 0,l,2,... 

where 5{q) is the phase shift experienced by particles scattered off the kink 

-1 3g„ 



-2 tan 



(1.15) 



(1.16) 



The allowed values merge into a continuum when L S> l/m. Thus, in the infinite 
volume limit, the ratio of the bound state energy to that of the lowest-energy scattering 
state is y^3/4 = 0. 



2. Quantum kink mass 

In this section we review existing approaches to calculating the kink mass. One can calcu- 
late the leading quantum correction to the classical kink mass in the weak-coupling limit 
A/m^ <C 1 analytically by summing up the zero-point energy contributions from all the 
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excitations (|lT^|lT4). This corresponds to one-loop level in perturbation theory, and 
gives §] 



Miioop = Mc\ + m 



(2.1) 



This result has recently been generalised to finite volume in Refs. | 

A non-perturbative alternative is to consider the response of the system to a 'twist' 
The mass of a topological defect is then defined as the difference between the free energy 
Ft-w in the topologically nontrivial sector containing the defect (in this case a kink), and 
the free energy Fq in the topologically trivial sector. 



1 
T 



(Ft 



tw 



1 1 Zo 
77; log— , 



(2.2) 



where T is the length of the system in the time direction. Unfortunately this is impossible 
to measure in Monte Carlo simulations because one only samples the partition function. 
We therefore resort to calculating derivatives of this mass with respect to a parameter g 
in the Lagrangian: 



1 

T 



(2.3) 



dg/tw \^9/o^ 

One then integrates this expression from a location where the mass of the defect is known 
exactly, such as the phase transition where the kink mass vanishes. In this paper we will 
take parameter g = m^, and then (^) becomes 



dm? 



L 
2" 



/tw 



(2.4) 



Unfortunately, linear error propagation is not very reliable for estimating the error in 
this integral; indeed, if we add more points this will tend to increase the error without 
any justifiable cause. One must also consider the error resulting from the quadrature. In 
practice, errors in the kink mass are under better control if one uses finite differences pO|| ; 
error propagation can then use the standard linearised results. We have then 



M{ml) - M{m\) 



1 



lU' 



e 2 



(2.5) 



mf ,0 



a formula which resembles the 'resampling' technique of Ferrenberg and Swendsen |]33| . 
Indeed, one way of checking that the measurement spacing is appropriate is to check that 
the measurements for the change of M from m\ to are the same when taken from above 
or below. 

An equivalent criterion, applicable to ( |2.4D , is to actually use the 'resampling' method 
to ensure that neighbouring values of {(fP') agree under the shift of m?. Unfortunately, this 
does not remove the problem with propagation of errors. 

One of the main drawbacks of this approach is the requirement to numerically integrate 
(or take finite differences) from, for example, the phase transition. Close to the phase 



transition, the errors in the expectation values in (2.4) will be quite large due to critical 
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slowing down. Worse, the overall error in the kink mass will increase with the length 
of the integration path. A local calculation of the mass of a topological defect will give 
better control of the errors and provide a check on the results obtained by the traditional 
approach. 

3. Correlation functions 

In this section we consider correlation functions in the topologically nontrivial sector. 
Specifically, we will work with two- and four-point functions of the scalar field. We use 
the two-point functions to study the kink mass, and the four-point functions to study the 
kink spectrum. These approaches are motivated by the spectral expansion of the field in 
the presence of a kink. 

3.1 Spectral expansion 

The correlation function of general operators Oi{t) localised in time has a spectral expan- 
sion 

oo 

Cijit2 - h) = {0,ih)0,it2)) = ^(0|Oi|a)(a|O,-|0)e-(*2-*i)^". (3.1) 

0=0 

where is |0) is the ground state of the relevant topological sector and \a) are a complete set 
of states with energies Ea- By a suitable choice of operators, we can determine individual 
terms in this expansion. By calculating the corresponding correlator in the presence of a 
kink, we can therefore determine its complete excitation spectrum, at least in principle. 
In practice, the energies can be determined either by a straightforward fit to 



or by using the Liischer- Wolff method |34, 35, 36|. In the latter approach, we consider the 
generalised eigenvalue problem 

= Xn{t,to)aj{to)p^^^ (3.2) 

where the eigenvalues A„ have the long-distance behaviour A„(t) = e~*'^" as t ^ co. The 
energies are then obtained with 



We show in Sections 3^ and 3^ that by considering operators Oi constructed from one 
or two field operators (j), we can measure not only the excitation spectrum in the rest mass 
of the kink, which corresponds to the linearised results in ( |1.12| - pn^ ), but also approximate 
wave functions of these states and the mass of the kink itself. 

3.2 Two-point function 

Let us first consider operators Ok = (l>{k) = f dx e^^^cf){x) for different momenta k, which 
are quantised because of the antiperiodic boundary conditions, 

k=Q^. (3.4) 
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Momentum conservation requires the momenta of the two operators to be equal, so for 
each momentum k we have the correlation function 



Cit2-ti;k) 



as calculated in our Monte Carlo simulation for the kink sector of the broken phase. This 
has a spectral expansion 



C{t2-h-k) = ^(0|(/>(ti;A;)|a)(a|(/)(t2;-A:)|0)e- 



it2-tx)Ec 



(3.6) 



Because of momentum conservation, all states \a) in this expansion must have overall 
momentum k. The lightest such state corresponds to the boosted kink. Other states in the 
expansion correspond to excited states of the kink and two-particle states consisting of a 
kink and a scalar particle. In the limit k <^ m, there should be a significant gap between 
the energies; the correlator (|3.5|) is therefore dominated by the boosted kink at long time 



separation t2 — ti. By analogy with the eigenspectrum given in Section 1.2, we expect 



C{t2 



ti; k) =(0|(/.(ti; fc)|kinkfe)(kinkfc|(/>(t2; -A;)|0)e-(*2-*i)^''"' + 
+ (O|0(ii; A;)|bs)(bs|</)(t2; -A;)|0)e-(*2-*i)^'^ + 

;(O|0(ii;A:)|a)(a|(/>(t2;-A:)|O)e-(*^-*^)^-'5' 



oo 
a=l 



(3.7) 



where is the energy of a boosted kink with momentum k; 



f and E^^^ are the 



energies for the bound state and a scattering state (a particle in the bulk) respectively. 
This spectrum has only one massive bound state, but we do not rely on this assumption 
in the calculations which follow. 

The physical interpretation of this is that the dominant contribution to the correlator 
comes from the (j) particle created at ti absorbed by the kink before re-emerging to be 
annihilated at t2- Not only must momentum be conserved by this process, but we also 
impose periodicity in the time direction: (j){x,t) = (t){x,T + t). As long as k <^ ni and 
(^2 — ti) ^ 1/m, the kink width is small relative to other relevant length scales, we can 
treat the kink as a point particle. Close to the phase transition, m approaches zero and the 
kink width diverges, and lower momenta k and longer time separations have to be used. 

Consider the action of a point particle of mass M in free space whose momentum 
increases by k at ti and decreases by k at t2 = ti +t. Periodic boundary conditions require 
that the initial and final x-components of the velocity must match. Then, 



S = M 



tJl + vi + {T-t] 



(3.8) 



noting that vq and vi are related by i <- 
relativistically, so fo, <C 1 and obtain 



(T — t). We assume that the kink moves non- 



(3.9) 
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and, by applying the stationary phase approximation we expect to see a dominant contri- 
bution to the correlation function, as measured on the lattice, that has the form 



C{t;k) 



Gs 



exp 



1 e t{T - 1) 

'2M T 



(3.10) 



In the infinite-time limit, T — > oo, this reduces to the non-relativistic kinetic energy of 
a particle with momentum A: [12|. In practice, the finite T effects are important, and 
therefore we keep the equation in this form. By fitting the asymptotic long-time behaviour 
of the correlator ( |3.5| ) to Equation ( p. 10 ), we can determine the kink mass M because the 
momentum k is fixed. This gives us a direct way of measuring the kink mass. 

3.3 Four-point functions and excitation spectrum 

To find the excitation spectrum in the rest frame of the 
kink, we need to use operators with zero overall momen- 
tum. Therefore we consider two-particle operators which 
correspond to creating two particles on a timeslice, sepa- 
rated by distance^ Ax. 

OAx{t) = ^(l){x,t)(l){x + Ax,t) (3.11) 



Figure 1: Schematic pic- 
ture of 'recoil' behaviour at 
a classical level. The kink 
is forced to match up due to 
periodic boundary conditions 
in the time direction (length 
T). It receives a kick from the 
scalar particle at ti which, to 
conserve momentum, it must 
give up at t2 — ti + t. 
0{L-^). (3.12) 




We can use the Liischer- Wolff method ( |3.3| ) to determine 
the energy spectrum {-Eq,}, which now includes only states 
in the rest frame of the kink and should therefore correspond 
directly to the negative parity states in ( |1.12| - pn^ ) . Whether 
a given state is a localised bound state of the kink or a 
free particle state can be determined by investigating its 
dependence on the volume L. In the infinite-volume limit, 
free particle states should behave as 



(2n + 1)2^2 



m 



In contrast, a localised bound state would have an exponentially small finite-size effect. 
When the volume L is smaller than the inverse mass 1/m we would expect the kink to be 
'squeezed' and the bound state to have lower energy. 

The other way of seeing the localisation of the first state is to try to reconstruct the 
wavefunctions from the generalised eigenvectors of (^]^). We can then compare our results 
with the eigenfunctions of the fluctuation field in the presence of a continuum kink (p..^), 
or its discrete equivalent discussed later. 

Using (3.1) to rewrite the generalised eigenvalue problem (|3.2|) 



E (0| Oax \l) e-^'* |0) = A. ^ (0| Oax \m) (m| Oa. |0) p^^l (3.13) 



(n) 



Ay 



^One could also work with Fourier-transformed fields with appropriate antiperiodic momenta. 
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we expect A„(t) = e " at long distance and so 

M OAy |0) p^^l = 5mn (3.14) 

Ay 

is a reasonable ansatz. Using our choice of two-particle operators this becomes 

^ H (t>{y)<P{y + Ay) |0) = 6mn- (3.15) 
y Aj/ 

We now split the field <p into the kink background (pk and fluctuations (p. We are in the 
broken phase; we assume in our expansion that the one-particle states are well-separated 
from the two-particle states. Then for the one-particle states 

Yl Yl ^"^1 ['^'^(y + Ay) + My - Ay)] 0(y) |0) = 6mn. (3.16) 
V 

We further assume that the field (f) can be decomposed into orthogonal modes - the wave- 
functions for each energy eigenstate. Noting that this means 

Y,^rn{x)(i)n{x) = 5mn (3.17) 

X 

we identify 

k{x) = YpaI IMx + Ax) + 0fc(x - Ax)] . (3.18) 

Ax 

This shows that, for one-particle states, the convolution of the generalised eigenvectors 
of the correlation matrices C{t), C(to) with the kink background gives the wavefunc- 
tion. A starting assumption for the form of the kink background in the weak-coupling limit 
is that it takes the standard shape (|1.3| ) but with the renormalised mass taking the place 
of m. This can be obtained from simulations with periodic boundary conditions. 

4. Numerical calculations and simulations 

For calculations on a lattice of size L x T, we make the standard choice of discretisation: 




where the summation over x runs over all sites, and /i is a unit vector in either direction 
on the lattice. For the remainder of this paper, we have set a = 1, and instead vary A and 

2 

m . 

Periodic boundary conditions are employed in the time direction at all times. The topo- 
logically nontrivial sector, where the kink is present, is simulated by imposing antiperiodic 
boundary conditions on the lattice in the spatial direction. 

Finite lattice spacing changes the dispersion relation by introducing a momentum 
cutoff. For a free particle, the dispersion relation is 

E(k) = ^-^sin^^+m^. (4.2) 
This modifies the 0{L^^) contribution in Equation ( p. 12 ). 
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4.1 Linearised calculations 



The linearised energy spectrum in Equations ( |1.12[ - |lT4 ) is calculated in continuum, and 
to be able to compare it with our Monte Carlo results, which are obtained on a discrete 
lattice, we need to understand what effect the discretisation itself has. We can see this by 
considering the Schrodinger problem for a kink in discrete space. 

A straightforward way of understanding the finite-size and discretisation effects is 
to solve the lattice Schrodinger equation for a kink background. Splitting into kink 
background (pk and fluctuations (j) a-s before, we find the classical discrete kink by the 
gradient descent method. 



(4.3) 



where antiperiodic boundary conditions are used and At is chosen to give good convergence. 
The eigenvectors and eigenvalues of (f) are then obtained by solving 



m A _ 
— 1 + -B 
2 4 



4>n{x) = E^(j)n{x) 



(4.4) 



where is the second-order lattice derivative with antiperiodic boundary conditions, and 
B is the diagonal matrix with entries B^^^^ = ^^(x)'^. 



0.2 



"a 0.15 




-B — B — B B &- 



0.05 



O D 0„ 



o.ooor 



0002 



Figure 2: Plot of the eigenvalues E"^ of the discretised Schrodinger equation (4.4) compared with 
the corresponding results for the continuum with a finite box, for A — 1/16, m? — 0.084. Plotted 
are the Goldstone mode (circles), the bound state (squares) and the first three scattering states 
(diamonds, triangles, crosses). There is a significant finite-size behaviour at this choice of m^, which 
will be worsened by radiative corrections. 
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The eigenvalues are compared with the continuum result ( 1.15 ) in Figure |^ for m 



2 



0.084, A = 1/16. These cannot be directly compared with the Monte Carlo results as the 
mass will pick up radiative corrections, but serve to give a qualitative description of the 
discretisation effects. 

4.2 Monte Carlo simulations 

Our Monte Carlo simulation uses a standard Metropolis algorithm, with acceptance rates 
at around 70%. A hybrid Monte Carlo algorithm was tried but not used for the results 
presented here, as it did not deliver a significant overall improvement in performance. Gen- 
erally, the space between measurements was kept longer than the autocorrelation time for 
the observable of interest. For the Liischer- Wolff eigenvectors, it is impossible to determine 
the autocorrelation time but it was assumed to be far longer than that for the two-point 
function results. 

4.2.1 Results for the kink mass 

As a benchmark, we first computed the kink mass using the conventional twist method 
(see Section ^) for many different values of m? along the line A = 1/16 in parameter space. 
The spatial average of cfP' was measured for antiperiodic and periodic spatial boundary 
conditions, so that the mass could be calculated from (p.5|). At least 2 x 10^ measurements 
were carried out for each value of m^, separated by 50 update sweeps. The integrated 
autocorrelation time for cjp' p7| , 

k2 a2 



^ 1 ^ {rnK+r}-{rn){rn+r) .^ 
'mt,</.2 „ 2^ /,4V _ / ,2\2 ' y^-^> 



was then estimated and used to find the number of independent measurements available 
for both the periodic and antiperiodic cases. The measurements were then binned appro- 
priately. These c/)^ measurements were then used to calculate 

/i = -lln(e-^(-i-?)E.0^)^^, (4.6) 

/, = _lln(e-^K-?)E.0A (4.7) 

in both the kink and trivial sectors. By analogy with the resampling technique of Ferrenberg 
and Swendsen, the measurement spacing has been made small enough that fi and /2 agree, 
within errors, in each sector. The final errors in the kink mass derivative were estimated 
by applying the bootstrap method. The kink mass is then obtained by using these data to 
calculate 

M{ml) - M{mi) = ^(/i,tw + /2,tw - /i,o - /2,o) (4.8) 

where 'tw' demotes the results of simulations with antiperiodic (twisted) boundary condi- 
tions and '0' those with periodic boundary conditions. The associated error is then 

A [M{ml) - M(m?)] = ^ [A/^^w + A/I,tw + A/'o + A/|,o 

+ (/i,tw - /2,tw)' + (/i,o - f2,of] . (4.9) 
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These expressions can then be used to find the kink mass 



N-l 

E 

n=0 



(4.10) 



Standard error propagation techniques can then be used to calculate the error in M{rn?). 

Our new method, outlined in Section can be used to calculate the mass at any 
point without depending on measurements at any other. We calculate the correlator (^ 
on the lattice, 



C{ti-t2]k) 




k 



(2n+ l)7r 
L 



and fit it to Ai, A2, M and E in the function 
C{t; k) = A, exp ' " 



+ A2 



^-Et^^-EiT-t) 



(4.11) 



(4.12) 



2M T 

where the first term corresponds to the contribution ( p. 10 ) from a moving kink and the 
second term to a free scalar particle. A fit to either a free scalar particle or ( 3.10| ) alone 
does not converge. The bound state of the kink cannot be determined by this fit, as it will 
be suppressed by a factor of 1/L. 

To exclude short-range behaviour, the fit is re- 
peated excluding more and more short-distance mea- 
surements until a 'plateau' for the fit parameters is 
reached. Our error estimates for the fit are obtained 
by performing an elaborate bootstrap: we resample our 
measurements of the correlator C(ti — t2',k), then re- 
peat the fit to obtain a series of estimates [^] . 

Our results, for varying m^, are shown in Figure 
^ Our results agree with [39|, in that the nonpertur- 
bative Monte Carlo result lies below the semiclassical 
mass result^. The two measurement methods are in 
agreement, from very close to the phase transition to 
beyond M ~ a = 1 where the treatment of the kink as 

a semiclassical particle might be expected to break down. The same number of measure- 
ments were made for each parameter choice in both methods. This means that the overall 
number of measurements needed for any given point is greatly reduced by studying the 
two-point function, as the errors are comparable in the two cases. 

In Table ^ we show the relative contributions of 'kink scattering' and 'bulk' behaviour 
to ( |4.1l| ) for a value of deep in the broken phase, and a sample plot of the corresponding 
data is given in Figure ^. 

4.2.2 Results for the particle spectrum in the presence of a kink 

We measured the operators ( p.ll ) for separations Ax = 0, 1 . . . ,L/2 — 1. An important 



Figure 3: Example Monte Carlo 
correlator data and fit for T = 48, 
= 0.1, A = 1/16 with k = 37r/48. 
Three points at each end are ex- 
cluded from the fit. 



requirement is that these operators are linearly independent; to ensure this we only use the 
^ There is an error in the equivalent calculation in |26| which was previously noted in |39|. 
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Figure 4: Comparison of tiie Icink mass measurements from a twist (Section g; shaded region with 
line at centre) and from the scattering method (Section crosses). The classical (L4) and first- 



order perturbative (2.1) results are also shown as solid and dotted lines respectively. Both types 
of measurement were carried out for L = 48; volumes L — 32 and L — 40 show finite-size effects 
that are similar to the magnitude of the bootstrapped errors. For the scattering measurements the 
errors increase when M > 1, but the technique is still valid. 



k = 


tt/L 








L 


Kink weight Ai 


Bulk weight A2 


E 


M 


32 


2.83875 ± 0.00095 


0.02293 ± 0.00081 


0.38382 ± 0.01332 


2.22514 ± 0.05441 


40 


3.02296 ± 0.00152 


0.01885 ± 0.00077 


0.39537 ± 0.05527 


2.29455 ±0.10858 


48 


3.13151 ± 0.00047 


0.01743 ± 0.00019 


0.41412 ± 0.00800 


2.19456 ± 0.03601 


k = 


3it/L 








32 


0.09402 ± 0.00041 


0.03040 ± 0.00052 


0.44331 ± 0.01196 


2.04267 ± 0.07687 


40 


0.14710 ± 0.00043 


0.02405 ± 0.00028 


0.42003 ± 0.00768 


2.13589 ± 0.07250 


48 


0.19189 ± 0.00044 


0.02016 ±0.00015 


0.41473 ± 0.00602 


2.18788 ± 0.04752 



Table 1: Contributions to the correlator (4.11) from 'kink scattering' and 'bulk' behaviour for 
TTi^ — 0.1 and A — 1/16 on a L x L lattice. The errors quoted here are from a bootstrapping method. 
At higher momenta the kink contribution becomes smaller, relative to the bulk contribution. How- 
ever, as the lattice size is increased the kink contribution increases. The kink mass remains within 
the error shown in Figure ^ for all measurements at fc = tt/L. 



first L/2 operators |35]. This is a similar imposition in position space to that of Liischer 
and Wolff in momentum space for their original paper [561. Our next step was to solve the 



generalised eigenvalue problem (3^) to give the energy spectrum. 

To find the maximum time separation t at which these measurements can be taken, 
we used the 'self-adjusting exponential fit' of Gockeler et al. [^], which compares sorting 
by largest eigenvalue at successive times {t, t + 1) to sorting by magnitude of the scalar 
product of eigenvectors. The two sorting procedures no longer agree at a time t where 
the data are too noisy to be reliable. Another check on the extent of the noise present is 
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the asymmetry of C/\x,Ay The noise in the data readily becomes apparent when fits or 
wavefunction reconstruction are attempted but these procedures serve as valuable checks. 

In order to verify the presence of a bound, localised state in the theory, simulations 
were carried out on lattices of spatial size L and time size 3L at various lattice sizes, to see 
if the squeezing of the kink and the bound state could be observed. In these simulations, 
A = 1/16 was used, with = 0.084 and m? = 0.504. These parameters have been chosen 
so that Xjrr? is slightly less than unity while 1/m is large enough to see the squeezing of 
the kink at smaller lattice sizes. This squeezing effect is seen in the ratio of the bound 
state energy to that of the first scattering state, plotted in Figure ||. 

The finite-size effects and the phase shift due to the kink can be seen when we attempt 
to extrapolate to the continuum limit (Figures ^ and 0). These should be contrasted with 
the Schrodinger results of Figure |^. 

For both plots, the lowest energy level is fit to a constant since it should not have any 
volume dependence. For the other excitations we can include the finite-size effects. From 



(3.12), this motivates a fit to c\ and C2 in 



2 , (2n + l 
m -|- ■ 



L2 



(4.13) 



for Figure |^. The value of re? is obtained from simulations with periodic boundary condi- 
tions. The data in Figure]^ are, however, too noisy to permit the same treatment. Instead 
we show the asymptotic behaviour given by the first two terms. 

4.2.3 Eigenvectors and wavefunctions 

In addition to the finite-size behaviour, further evidence for the bound state is given by the 
shape of the eigenvectors. In Figures |^ and ^, the eigenvectors of the correlation matrix 




Figure 5: Plot of the ratio of the first to the second energy eigenvalue against lattice size, for 
A = 1/16, rtT? = 0.084. The classical continuum result mhs/iri^ is shown (dotted line), as well as 
the estimated continuum result obtained by the intercepts of Figure ^ (solid line) . The results here 
are at t = 3. 
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0.8 



0.4 - 



0.2 ■ 



Figure 6: Plot of the first five energy levels given by the diagonalisation of (3J_) for the operators 
( 3.11 ), with A ~ 1/16, TV? = 0.084. This was done at relatively short distance (t = 2). For the 
scattering states, the quality of the data here has allowed a constrained fit with two orders of 
discretisation effects ( 4.13| ). 




0.0001 0.00015 



0.0002 



0.00025 



Figure 7: As Figure |^ but with A = 1/16, n? = 0.504. Here the lines are not fits, but show the 
asymptotic free-particle behaviour E'^ = + (2n + l)^7r^/L-^. 



are given for the first four energy levels at L = 64 (the error bars are estimates based on a 
bootstrap of the eigenvectors for resampled sets of measurements) . 

Following the method given in Section these eigenvectors can be used to reconstruct 
the original wavefunctions. First, we estimate the renormalised mass by fitting the zero- 
mode energies in the periodic sector to a constant. This renormalised mass is then used in 



the classical gradient flow (4.3) method to estimate the background kink. 

Approximate errors are recalculated using the bootstrap method; interestingly, the 
errors are much smaller than those for the corresponding eigenvector. The results of the 
reconstruction are shown in Figures M and 0. In Figure the kink is slightly wider than 
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the lattice spatial size, and hence the bound state is 'squeezed'. 



5. Conclusions 



We have shown that correlation functions in the presence of a kink provide detailed informa- 
tion about the kink itself: its mass and its excitation spectrum, along with the approximate 
wavefunctions. One can measure these correlation functions using standard Monte Carlo 
techniques, thereby calculating its properties in a fully non-perturbative way in quantum 
field theory. 

Quantum kink masses have been calculated previously in several different ways: at 
one-loop level in perturbation theory, and non-perturbatively using creation operators or 
by calculating the excess free energy of the kink across the two topological sectors. Our 
method, in which the mass is obtained from the kinetic energy of a kink with a known 
momentum, is closer in spirit to the way other masses are calculated in Monte Carlo 
simulations, and potential errors should be better controlled. 

Calculations of kink excitations have so far been restricted to the linear level, and it 
is quite likely that in many cases the interactions will modify even qualitative features of 
the spectrum, such as the number of bound states. It is therefore important to have a way 
of measuring the spectrum non-perturbatively. As a side product, we obtain approximate 
quantum wave functions of the energy eigenstates, which are only valid at weak coupling 
but which nevertheless help us to identify the states. 

We demonstrated our methods by carrying our Monte Carlo simulations at weak cou- 
pling in the 1 -|- 1-dimensional Xcj)'^ theory. The results show that the mass and the exci- 
tations can be determined in practice, with relatively small errors. It will be interesting 
to extend these simulations to the critical region, and to see explicitly how the behaviour 
changes when critical phenomena become important and expectations based on linear the- 
ory become invalid. 

It should be straightforward to generalise our methods to other topological defects and 
more realistic theories, by considering correlators of suitable operators. Particularly inter- 
esting applications of this would be dyons as excitations of 't Hooft-Polyakov monopoles, 
and localisation of degrees of freedom on cosmic strings, domain walls or solitonic branes 
in braneworld models. Comparing simulation results with known exact results would also 
provide a useful test of lattice formulations of super symmetric field theories. 
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Figure 8: Plot of four lowest-lying eigenvectors for L = 64, = 0.084. 
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Figure 9: Plot of four lowest-lying eigenvectors for L = 64, = 0.504. 
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Figure 10: Plot of four lowest-lying wavefunctions for L = 64, m = 0.084. 
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Figure 11: Plot of four lowest-lying wavefunctions for L = 64, = 0.504. 
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